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A regularization method is presented to deduce dynamic 
correlation functions from exact diagonalization calculations. 
It is applied to dimer-dimer correlation functions in quantum 
spin chains relevant for the description of spin-Peierls systems. 
Exact results for the XY model are presented. The analysis 
draws into doubt that the dimer-dimer correlation functions 
show the same scale invariance as spin-spin correlation func- 
tions. The results are applied to describe the quasi-elastic 
scattering in CuGcOa. 



I. INTRODUCTION 

Dynamical correlations characterize the spectral prop- 
erties of physical systems. They are accessible by a multi- 
tude of experimental setups. Even static properties such 
as correlation lengths extracted from energy integrating 
methods such as X-ray scattering also require the theoret- 
ical knowledge of the dynamics of the correlations to al- 
low for comparison. The access to dynamical correlation 
functions for physically relevant systems is usually diffi- 
cult even in exactly solvable models.uTJ Dynamical spin- 
spin correlation functions in..Heisenberg chains have b£fP 
widely studied numericallyoB as well as analytically.Brj 

Usually the focus lies on the imaginary part of the cor- 
relation functions. The real and the imaginary parts can 
be Kramers-Kronig transformed into each other and thus 
hold the same information. The information that can be 
extracted from finite systems accessible by exact diago- 
nalization (ED) concerning the thermodynamic limit is 
limited. The accuracy of the results is different for differ- 
ent regions of the spectrum. In the case of finite systems 
it proves thus useful to actually calculate both quantities 
to extrapolate to the thermodynamic limit. The empha- 
sis of this article lies on the presentation of a method to 
accurately extract the real part of the correlation func- 
tions. 

The random-phase-appraximation (RPA) approach 
applied by Cross and Fisherfl to describe the spin Peyk 
transition has been successfully applied to CuGeO; 
To compare with neutron or X-ray scattering data the 
phononic dynamic structure factor has to be determined. 



5(q,w) 



l E„hnA,(q,a;) 

7T 1 — exp(— [3uj) 



(1) 



mal coordinate propagator D u (q,oj) has been calculated 
including the spin-phonon coupling term in RPA and de- 
pends crucially on the dimer-dimer correlation function 

x(q,w). 
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The spin-phonon coupling constants are denoted g v > q 
and the unperturbed normal coordinate propagator is 
D„ (q, oj) . Especially the real part of x(q, u>) is of rel- 
evance since it determines dqminantly the poles of the 
phononic correlation functionli!J'EirE3 This dependence is 
the initial motivation for an accurate determination of 
the dimer-dimer correlation function as presented here. 

The RPA-calculations predict for CuGeC>3 the appear- 
ance of spectral weight in the center of the phonemic spec- 
trum as a precursor of the phase transition.t2lll3 I use 
the numerical results to discuss the temperature de-pep. 
dent intensity of the predicted precursor in neutronuJ't3 
and X-rayEZI scattering experiments. Also, the resulting 
description xtf the phonon hardening is consistent with 
experiment £3 

The relevant spin-phonon coupled Hamiltonian for 
spin-Peierls systems consists of three parts H = H s + 
Hp + H sp . Of relevance here is the Heisenberg spin-chain 
Hamiltonian 
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with the superexchange integrals J and J2 between 
nearest-neighbor (NN) and next-nearest-neighbor (NNN) 
magnetic ions, respectively, and spin 1/2 operators Si at 
magnetic ion site 1 in the three-dimensional lattice, z is 
a unit vector along the spin-chain direction. For a de- 
tailed discussion of the harmonic phonon part H p and 
the spin-phonon coupling term H sp please refer to Ref. 

n 

Defining the Fourier transform of the dimer operators 



(4) 



v labels the relevant phonon modes. The retarded nor- 



the dimer-dimer correlation function can be written as 
X(>k,iun) = ~ J ' dr e**»T (r q (r)Yl q (0)) (5) 



1 



with Matsubara frequencies ui n = 27rn//3, inverse tem- 
perature (3 — l/(k-aT), and number of unit cells N . It 
depends only on momenta q z along the spin chains. It 
has been calculated in the analytically continuted form- 
where iu) n — ► u> + it with e — ► 0, by Cross and Fisheru 
with bosonization techniques as 
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with the spin- wave velocity v s , lattice constant c along 
the magnetic chains, and the functions 
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The result has the general form of spin correlation func- 
tions obtained from conformal field theoryBEl I set % = 1 
in this paper. 

The prefactor Xq^bT j J) is assumed constant in field 
theory but has been shown by Raupach et al. using den- 
sity matrix renormalization group (DMRG) studies to 
be temperature dependent in the static case and for 
q z = 7r/c,t3 where c is the lattice constant along the 
magnetic chains. This wave vector describes the lattice 
modulation in the spin-Peierls phase and is thus the rele- 
vant one for this paper. I will now turn to the question of 
how far the dynamical properties of the correlation func- 
tion are correctly described by the field-theoretical result 
Eq. (D. 



II. REGULARIZATION 

The correlation function can be calculated after the 
diagonalization of the spin Hamiltonian in the spectral 
representation since eigenfunctions \n) and eigenvalues 
E n are known. All numerical results in this paper are 
obtained using periodic boundary conditions. Defining 
the matrix elements 



V nm = (n \Y qz | m) 



and the Boltzmann factor 
1 
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where Z is the partition function, one can write 

(P) | V nm | (lo + E n — E„ 
Re X (q z ,Lo) = lim ^ _ F ,2^,2 
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Since real and imaginary part of the correlation function 
are Kramers-Kronig related via 
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it suffices in principle to know one of them, or, in other 
words, both parts contain the full information on the 
correlation function. Since the diagonalization of the 
Hamiltonian usually is limited to finite system sizes the 
information that can be extracted concerning the ther- 
modynamic limit is limited. In that case it proves use- 
ful to actually calculate both quantities. Note that for 
IuY\{q z , ui) to have a cut off, i.e., the correlation function 
has a finite spectrum, Re x(9zj w) must change sign as a 
function of lo. 

The spectra of finite systems are discrete and thus the 
correlation functions Eqs. (10) and ([ll]) contain a series 
of poles as a function of uo at frequencies uo nm = E rn — E n . 
In order to extract information on the correlation func- 
tion in the thermodynamic limit, i.e., N — > oo, one has 
to distinguish between well defined quasi-particle spectra 
and those where interaction induces transitions between 
many states. In the XY case excitations of the non- 
interacting spin-less fermions form a well defined con- 
tinuum in the N — > oo limittj as well as the singlet exck 
tations out of the ground state in the Heisenberg chain.El 



A. Imaginary part 

The averaging via binning is obtained by replacing the 
delta function in Eq. (|ll|) by a step function with an area 
of unity. 
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9(x) is the Heaviside function. For systems with well de- 
fined quasi-particle excitation spectra with spectral lines 
at frequencies Coj as the XY-model or Heisenberg chains 
at T = 0, the appropriate choice is such that the interval 
boundaries lie in the middle between those spectral lines: 
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Introducing the density of states rj~ 1 with resp 
Bcthe-ansatz quantum numbers it can be shownuu 
that for systems with well defined quasi-particle excita- 
tion spectra Eq. (|ll]) can be represented scaled with the 
system size as 

uij=E n -E m 

Imxfe>^) = - lim —T}[Qj] V fnm((3) \V nm \ 2 . 

I, )fYl / ' * 



(16) 



2 



The sum covers only values of n and m such that ujj — 
E n — E m , L is the number of magnetic ions on a chain. 
Approximating the derivative in the density of states for 
a finite system as 



1 u>i — £j,j 
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Eq. (Ill) becomes 



£j i= E n -E„ 



E 



2n/L 



-2tt f nm (f3) \V m 

UJj+l — 



(17) 



(18) 



Given the definitions Eqs. (I14J) and (|15|) the expres- 
sions ( [i"3| ) and (jl^) are identical, except for the latter be- 
ing only defined at discrete frequencies. Both approaches 

■yco] 

( pL5|) in one case and of 
in the other. For larger 



suffer the disadvantage of ambiguit y c oncerning the defi- 
nition of the interval Eqs. ( |l4|) a nd (| 

fI3 



the difference, .quotient Eq 



system sizesJ3 ! E3 where the density of spectral lines ujj is 
much higher, this ambiguity becomes irrelevant. 

For Heisenberg chains at finite temperatures the defi- 
nition of the density of states rf 1 with respect to Bethe- 
ansatz quantum numbers is inappropriate since a large 
number of spectral lines are present (see Sec. IV). The 
binning procedure Eq. ( |l3| ) has to be applied. At higher 
temperatures the interval bojindaries 
ally are chosen equidistantly.Bca 



B. Real part 
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As discussed in the subsequent sections it turns out 
that for Auj > 0.2 J Hex(Qx,&j) is a smooth function 
of Aw. For the Heisenberg chains discussed Re x((Zz, &j) 
is constant for 0.06J < Auj < 0.13J. At intermediate 
temperatures 0.1J < k^T < 0.9J the method yields in- 
accurate results for Auj < 0.6 J since the data points are 
located on shoulders of nearby divergences. At higher 
temperatures the finite size effects become negligibleJEI at 
lower temperatures the limiting value is given by the ac- 
curacy with which the ujj are determined. I show that a 
choice of Auj = 0.1J yields reliable results in all param- 
eter regimes. 

For higher frequencies the results for the real part of 
the correlation functions are free of finite size effects. 



III. XY MODEL 

I demonstrate the regularization method and its ap- 
plicability for an exactly solvable case, the XY model. 
Neglecting the exchange between the Sf components of 
the spin system the magnetic Hamiltonian reduces to 



(21) 



I have also neglected the next nearest neighbor coupling 
Ji- The dimer operator becomes 



(22) 



As can be seen from the denominator of Eq. ([10|) the 
real part of x is divergent for frequencies uj — > E m — E n . 
For arbitrarily small but finite e the term which leads 
to the divergence at uj ~ oj nm = E m — E n vanished for 
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It turns out that for systems with well defined quasi- 
particle excitation spectra as the XY-model or Heisen- 
berg chains at T = and for an accordingly chosen set 
of frequencies {ujj} — {u) n m} this allows a very accurate 
determination of the real part of the correlation function 
from Eq. (jTo)). 

In the case of Heisenberg chains at finite temperatures 
additional spectral lines stemming from transitions out of 
various thermally activated states are close to each other. 
This results in lack of accuracy since the frequencies uij 
chosen at T — now lie uncontrolled on the shoulders of 
those additional spectral lines. This problem can be over- 
come by suppressing all contributions to the correlation 
function value in an interval E n — E m G [ujj ± Auj] . 
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A. Exact results 

The spin operators in this model can be transformed to 
non-interactingpSpin-less fermions via a Jordan- Wigner 
transformation.BHI The correlation function-can then be 
determined in the thermodynamic limit J13 ! 

™ 9.T, 2—i uj+E,, — Ei. , _ -I- if 
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The energy dispersion is given by Ek z = J cos(k z c), 
are Fermi distribution functions, the sum covers the first 
Brillouin zone, L is the number of magnetic ions in a 
chain. In general for L > 10 4 the result is independent 
of L for all practical purposes. 



The imaginary pan^pf Eq. (23) can be given for q z 
tt/c in closed form asc3 



Im X XY ( 7r / c ' w ) : 
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In the limit T — > the real part attains a simple form, 
too. Defining x = uj/(2J) one has for \x\ < 1 



3 




FIG. 1. Dimer-dimer correlation function in the XY case 
at T = 0. (a) imaginary part, (b) real part. Thermodynamic 
limit results (full lines) are obtained via Eqs. (|2j| ) through 
. Step functions are binned as given in Eq. (UjJ. Symbols 
are from Eq. ( |Ti| ) (imaginary part) and ( pjj ) (real part) at 
frequencies uij. Finite size effects of the real part are largest 
at uj — 0. Dash-dotted line: exact diagonalization from Eq. 



FIG. 2. (a) Full line: imaginary part of the correlation 
function Eq. (124) at feT = 0.3J. Step functions are binned 
as given in Equity. Symbols are from Eq . (jig) , (b) Full line: 
real part of the correlation function Eq. (ptyat ksT = 0.3J. 
Dash-dotted line: exact diagonalization from Eq. Jit)), Full 
symbols are values at the spectral lines ujj £ Wxy present at 
T = 0, open symbols are additional spectral lines VV^y \ Wxy 
at Snite T. 



and for Ixl > 1 
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(25) 



Rexj^(7r/c,w) = ^ X arcsin(l/x) . (26) 

The imaginary part of the correlation function Eq. fl23| ) 
is depicted for different temperatures by the solid lines 
in Figs. |l| and || (a), the real part in Figs. |l| and || (b). 



B. Comparing finite and infinite systems 

In finite systems with periodic boundary conditions 
the allowed wave vectors are k z (j) = 2ttj/L for L/2 
odd and k z {j)rfF n{2j + 1)/L for L/2 even with j S 
{0, . . . ,L - 1)£5 For q z = ir/c at T = the only spec- 
tral components contributing to the dimer-dimer cor- 
relation function Eq. (|2^) are fermion-hole excitations 
with ujj = —2E k ^jj forming a well defined continuum 
in the thermodynamic limit and wi±Jb-excitations only 
within the half filled fermionic band.BEj In the example 
of L = 14 the values to consider are u)j £ Wxy > where 

w xy = ^{0-445, 1.247, 1.802, 2}. 

In Fig. [y I show different approaches to deduce the 
imaginary and real part of the dimer-dimer correlation 
function from finite systems compared with the exact re- 
sult in the thermodynamic limit [L — > oo, full lines) at 



T = 0. The step functions in Fig. [j] (a) resulting from the 
binning approach Eq. ( |l3| ) approximate the exact curve 
quite satisfactorily, albeit with limited freq uency reso- 
lution. The scaling approach from Eq. (18) also yields 
excellent agreement as shown by the symbols. 

The dash-dotted line in Fig. [j] (b) is the result of the 
real part from Eq. (10) via exact diagonalization for 14 
sites. Above the fermionic spectrum, i.e., u> > 2J, the 
finite system gives a very accurate result. This is equiv- 
alent to the good agreement of the real time correlation 
functions from ED and L — > oo results for small time 
scales^ The symbols show the results using Eq. ( |l0| ) by 

only regarding frequencies W X y an ^ £*> = 0. The points 
give a satisfactory approximation of the exact curve in 
the thermodynamic limit. The discrepancies due to fi- 
nite size effects are largest at to = 0. 

At finite temperatures additional spectral lines are 
present, namely for all wave vectors k z (J) — irj/L, al- 
beit with temperature dependent weight. The additional 
excitations are fermion-hole like out of thermaly excited 
superpositions of states with a fermionic vacancy and 
an additional fermion |$±) = \/2 _1 (|L/2 - 1) ± \L/2 + 
1)). For L = 14 on has Qj e W^y with WgV = 
J{0.445, 0.868, 1.247, 1.564, 1.802, 1.950,2}. This cir- 
cumstance can be used to enhance the energy resolu- 
tion of the calculations at finite temperatures. Since the 
spectral weight of the additional spectral lines is small 
at small temperatures, an application to the imaginary 
part will only yield reasonable results for T > 0.3J. Fig- 
ure U (a) shows the excellent agreement with the exact 
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curve (full line) of both the binned step functions from 
Eq. (|l3|) and the symbols of the scaling approach Eq. ( |is| ) 
using u {0}- Note that the application of the scal- 

ing approach is only possible here since I can take into 
account all spectral lines, uj = needs to be included 
in the spectrum because of the degeneracy of |$+) and 

l<M ' 

The calculation with the extended spectrum VVxy can 
be applied to the real part in principle at any finite tem- 
perature, albeit with some small scattering of the data 
points at k B T < 0.3J as shown for the Heisenberg cases 
in the subsequent sections. The continuous line in Fig. 
|2| (b ) shows the real part of the correlation function Eq. 
(||) at k B T = 0.3 J for the XY model. The dash-dotted 
line is the result from Eq. ([h]) via exact diagonalization 
for 14 sites showing the spectral lines tlij as divergences. 
For us > 2J finite size effects are absent. The full symbols 
show the results using Eq. ( |To| ) by only regarding frequen- 
cies as present at T = 0, open symbols are added at 
finite temperatures. The points give an excellent approx- 
imation of the exact curve in the thermodynamic limit. 



IV. HEISENBERG MODEL 

I now determine the dimer-dimer correlation function 
for Heisenberg chains with the Hamiltonian given in Eq. 
(||) in the case of a — J 2 / J = 0. 

A. Zero temperature 



At T — the spectral lines u>j € Wg^ at q z = ir/c can 
be identified as the magnetization conserving (singlet) 
excitations out of the ground state determined via the 
Bethe ansatz, which form- a well defined continuum in 
the thermodynamic limit .Era Equation (18) can thus be 
applied and Eq. ( |l0| ) can be regularized as discussed in 
Sec. 0. 



For a 16 site chain the four spectral lines at Wg 16 '' = 
J{0.446, 1.669, 2.588, 3.083} lead to the imaginary part 
of the dimer-dimer correlation function as shown in Fig. 
H (a). The bins are obtained using Eq. ( fl3"| ) and the 
symbols are given using Eq. The_fLeld theoretical 

prediction Eq. (H) reduces for T — > toBTj 



Imxgkir/c ) w) = -0(A-M) 



(27) 



From the fit (dashed line in Fig. (3| (a)) I conclude C~l. 
The correspondece of the fit is not too good and suggests 
the presence of correction terms. The cut off parameter 
A rs 7r J is determined via the real part (see below) and 
indicates that the speqtr.al. weight lies mostly within the 
two-spinon continuum P 6 !! 2 *! 

The real part of the dimer-dimer correlation function 
is shown in Fig. |3| (b). For lo > ttJ the numerical results 
are in the thermodynamic limit showing the cut off to 




FIG. 3. Imaginary part (a) and real part (b) of the 
dimer-dimer correlation function in the unfrustrated Heisen- 
berg chain at T = 0. The imaginary part for Enite systems is 
binned (Eq. ( |J 3j) ), each bin holds one spectral line, symbols 
are from Eq. Jl8j) . The dashed line is the field theoretical re- 
sult (cf. Eq. fciW). The symbols for the real part as well as the 
dash-dotted line are obtained using Eq. (lid), for us > 3. 14 J 
the numerical result is in the thermodynamic limit (full line). 



be A = it J. For u> = finite size effects are significant 
since one expects from field theory and by analogy to the 
XY model Re x (0) (tt/c, uj = 0) -> -co. For u ^ the 
field theoretical result obtained by the Kramers-Kronig 
transform of Eq. (|27|) is 
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with a non-continuous jump to Re Xcp(tt/c, 0) = — oo. In 
contrast to that for ui < J the numerical data in Fig. || 

(b) indicate Re x (0) (tt/c,oj) < Rex^Tr/c, ui) < 0. This 
suggests significant corrections to the field-theoretical re- 
sults Eq. (|!) and @. 



> 



(28) 



B. Finite temperatures 

At finite temperatures additional spectral lines attain 
significant weight as can be seen by the divergences 
in the real part of the correlation function Eq. (|Io| ) 
as plotted for k B T = 0.1J in Fig. | (b) (dash-dotted 
line with L=16). The most dominant spectral lines are 
W^ 16) = J{0.446, 1.125, 1.669, 2.174, 2.588, 2.870, 3.083}. 

The imaginary part of the dimer-dimer correlation 
function is determined via the same binning procedure 
as for T = and shown in Fig. | (a). At k B T = 0.1J 



the weight of the additional spectral lines >Vq 16 ^ \ Wq 
is too small to improve the frequency resolution as has 
been done in the XY case for k B T = 0.3J. The scaling 
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FIG. 4. Imaginary part (a) and real part (b) of the 
dimer-dimer correlation function in the unfrustrated Heisen- 
berg chain at IcbT — 0.1J. The imaginary part for Enite 
systems is binned as in Fig ^. The real part for uj < 3. 14 J 
is obtained using Auj = 0.1J and is given by full symbols at 
frequencies as for T = and by open symbols at additional 
spectral lines appearing for T > 0. For u > 3.14J the numer- 
ical result is in the thermodynamic limit (thick full line). The 
full lines are the Cross and Fisher result Eq. the dashed 
lines are rescaled according to Eq. ( p^ ). 

regularization Eq. ( p"8| ) is inappropriate at finite temper- 
atures in the Heisenberg model as discussed in Sec. O. 

The real part of the correlation function has to be reg- 
ularized according to Eq. (|2^). The resulting values for 
Auj = 0.1J are shown by the symbols in Fig. |4| (b). Full 

symbols are evaluated at the Uj € Wq 16 \ open symbols 
are determined at additional dominant spectral lines and 
are less accurate. The finite size effects are of the order of 
10% at u> — as can be seen by comparison with DMRG 
results in the thermodynamic limitEj shown by the cross 
in Fig. ^ (b). They are negligible for uj > nJ. 

The full lines in Fig. [| (a) and (b) show the field- 
theoretical result Xcf(1z,^) as given by Eq. @. The 
functions are scaled by Xo(¥b\T '/ 'J) to match the uj = 
values obtained by DMRG.E3 Both real and imaginary 
part fit the numerical data better when rescaling the 
frequency dependence as x CF (q 2 ,w — > 0.54w) (dashed 
lines in Fig. ^ (a) and (b)). The temperature depen- 
dence of the discrepancy between field-theory (dashed 
line) and numerical result for the real part (full line) at 
uj = uji = 0.446 J is shown in Fig. ||. 

As demonstrated in Fig. the numerical result for 
Rex(-7r/c, u>i) with u)\ = 0.446J is practically indepen- 
dent of Aw for all temperatures. I use this circumstance 
to determine the scaling function such that 



R-eXcF^/C'ST^i) =Re%(7r/c,wi). 



(29) 
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FIG. 5. Temperature dependence of Rex(7r/c, (jJi) for 
uji = 0.446J: the result from Eq. ( |2c| ) is practically inde- 
pendent of Auj, the lines for Auj = 0.02 J and Auj = 0.1J fall 
on top of each other. The dashed line is the held-theoretical 
result ReXcF('/zi a; i) from Eq. m). The inset shows the scal- 
ing function defined in Eq. (bg) for different chain lengths. 
The frequency corresponding to L — 14 is uj\ — 0.516J, for 
L = 12 it is uj! = 0.61J. 



slope of the fit function as can be seen in Fig. [| (b) . The 
resulting gj> is shown in the inset of Fig. 0. The results 
for different chain lengths show that L = 16 yields an 
upper bound for gx, the finite size effects decrease with 
decreasing temperature. 

The following arguments show the validity of the real- 
part results even for finite regularization parameters Aoj: 

(i) All symbols shown in Fig. ^ (b) follow a smooth 
line as expected by comparing with field-theoretical re- 
sults and the XY model. The real part changes sign as a 
function of frequency as expected. 

(ii) The dependence of Re x^ (i"/c, ojj) on Auj is small. 
This is shown for ujj = oj\ = 0.446J in Fig. ^, where the 
the values of the correlation function for Auj — 0.02 J and 
Aoj = 0.1J fall on top of each other for all temperatures. 

(iii) The scaling function gx determined by 
Rex(7r/c,wi) also yields a better field-theoretical fit of 
the imaginary part of the correlation function. 



V. FRUSTRATED HEISENBERG CHAIN 

For J2 7^ the Heisenberg model is not integrable any 
more and thus the classification of the excitation spec- 
trum at T — is not possible from the Bethe ansatz. I 
discuss here a value of a — J2JJ1J 0.24 which is among 
the proposed ones far CuGeOsErEH and is close the quan- 
tum critical pointed making the field-theoretical result 
given in Eq. (g) eligible for comparison. 



A. Zero temperature 

The dash-dotted line in Fig. || (b) shows that at 
T = for a 16 site chain there are four spectral lines 



This procedure shows little finite size effects since the 
finite size scaling displaces the points at u>t along the 



(16) _ 



J{0.233, 1.174, 1.963, 2.403} which, by 
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a=0.24, k B T=0 
0.9 of' 




ca [J] 




FIG. 6. Imaginary part (a) and real part (b) of the 
dimer-dimer correlation function in the frustrated Heisenberg 
chain at T = with a — 0.24. The imaginary part for fi- 
nite systems is binned (Eq. (JI^)), each bin hold one spec- 
tral line, symbols are from Eq. (It ). The dashed line is the 
field-theoretical result (cf. Eq. (21)). The symbols for the 
real par t as well as the dash-dotted line are obtained using 
Eq. (Ud), for u > 2.5J the numerical result is in the thermo- 
dynamic limit (full line). 



FIG. 7. Real part of the dimer-dimer correlation func- 
tion for a = 0.24 at k B T = 0.1J (b) and k B T = 0.3J (a). 
Dash-dotted lines: results applying Eq. Symbols: val- 

ues obtained via Eq. ( |2C[ ) with Aui — 0.1J for different system 
sizes. Full lines are the Beld-theoretical results from Eq. (Q), 
the dashed line is the scaled result, see Eq. 



B. Finite temperatures 



analogy to the unfrustrated case, may also be identi- 
fied as the singlet excitations out of the ground state. 
It is thus reasonable to suppose them to form a well de- 
fined continuum in the thermodynamic limit and thus 
Eqs. (|l|) and @ can be applied (Au = 0). Note that 



W\ _ |v*,(i)| _ 



\W [ >\ = [£-Lmod(4)]/4 for the chain lengths 
investigated. 

The imaginary part of the dimer-dimer correlation 
function is shown in Fig. |^ (a). The bins are obtained 
using Eq. 



and the symbols are given using Eq. (18) 
The fit to the field theoretical prediction Eq. 
line in Fig. || (a)) leads to conclude C = 0.9 
parameter A 



|27| ) (dashed 
The cut off 

2.5 J determined via the real part (see be- 
low) indicates that the spectral weight lies mostly withm 
a bounded continuum just as in the unfrustrated case.Ea 
The real part of the dimer-dimer correlation function 
is shown in Fig. ^ (b). For uo > 2.5 J the numerical results 
are in the thermodynamic limit suggesting the cut off to 
be A = 2.5 J. For u> = finite size effects are significant 
since one expects from field theory and by analogy to the 
XY model Re x (0) (tt/c,uj = 0) -> -oo. 

The agreement with the field-theoretical results Eqs. 
( p7j ) and ( p8| ) is better than in the unfrustrated case 
but the analogy to the XY model again suggests 
Re x < '°H 7r / c, 0.233 J) < 0. It cannot be entirely excluded 
though that the discrepancy is due to finite size effects. 



At finite temperatures additional spectral lines attain 
significant weight as can be seen by the di verg ences in 
the real part of the correlation function Eq. (|0|) as plot- 
ted for k B T = 0.1J in Fig. fj] (b) and for k B T = 0.3J in 
Fig. ^ (a) (dash-dotted lines with L=16). The most dom- 
inant and thus relevant spectral lines are ojj € VV*^ = 
J{0.233, 0.720, 1.174, 1.616, 1.963, 2.232, 2.403}. Again, 
the analogy to the unfrustrated Heisenberg chain is ob- 
vious since IVV^I = IVV^I for the chain lengths inves- 
tigated. 

The imaginary part of the dimer-dimer correlation 
function is determined for k B T — 0.1J via the same 
binning procedure as for T = and is shown in Fig. 
|8| (b). For k B T — 0.3 J the set of frequencies given by 
U {0} enhancing the resolution as in the XY case. 
The resulting step functions from the binning procedure 
Eq. ( |l3| ) are shown in Fig. || (a). In both cases the scal- 
ing regularization Eq. ([l8]) is inappropriate as discussed 
in Sec. 

The symbols in Fig. [7| (a) and (b) show the numerical 
results of the real part of the dimer-dimer correlation 
function applying Eq. @ for k B T = 0.3J and k B T = 
0.1J, respectively, with Au — 0.1J. The full symbols are 
at spectral lines and at ui — 0, the open symbols 

are at W $ \ W $ . 

For ui > 2.5J the numerical results in Fig. |t] are 
in the thermodynamic limit. The multitude of spec- 
tral lines appearing at finite temperatures which lead 
to the requirements on Aw can clearly be seen by the 
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to [J] 

FIG. 8. Imaginary part for k B T = 0.1J and k B T = 0.3J 
in the frustrated Heisenberg chain. Binned curves from Eq. 
(pl[), dashed lines from field theory Eq. (0), the dashed line 
is rescaled as in Fig. [?]. 

dash-dotted lines in Fig. [7] (a) and (b). To obtain the 
smooth results of the regularized data points a choice of 
0.06 J < Auj < 0.15 J is necessary at intermediate tem- 
peratures 0.1J < k B T < 0.9J. 

The temperature dependent effect of the choice of Auj 
is shown in Fig. |9| for three different frequencies. At 
LU3 = 2.8J the result is equal to the thermodynamic limit. 
At intermediate frequencies, as shown in the example of 
lo 2 = 1.96J and uj\ = 0.23J in Fig. ^ (b) and (c), respec- 
tively, values of Aw < 0.06 J and Alu > 0.15 J deviate 
from the best choice of Alu < 0.1J. For other tempera- 
tures, i.e., 0.1J > k B T and k B T > 0.9J, one may choose 
Alu as small as 0.01J. The influence of finite size effects 
can be estimated by comparing the to — .-values with 
DMRG results shown by crosses in Fig. RJIj Together 
with the above discussed small dependence of the values 
on Aw this is strong evidence for the accuracy of the 
results. 

The field-theoretical result in Eq. (|J) is expected to 
well describe the low energy physics. The full lines in 
Fig. show the corresponding curves. The functions are 
scaled by Xo&bT / J) to match the u> — values obtained 
by DMRG.E3 While at k B T = 0.1J both the imaginary 
and the real part of the dimer-dimer correlation function 
are well described for to < J as shown in Figs. |s| and 
|?] (b), at k B T = 0.3 J the discrepancy to the numerical 
result is obvious as shown in Figs. || and fj] (a). 

As in the unfrustrated case I attempt to improve the 
field-theoretical result by introducing the scaling func- 
tion Qt defined in Eq. (|29|). The relevant frequency in 
the frustrated case is u>i = 0.23J for L = 16. The tem- 
perature dependence of x(7r/c, wi) is depicted by the full 
line in Fig. |^ (c) in comparison to field theory (short- 
dashed line). The scaling function gx is shown in the 
inset of Fig. (c). The dashed line in Fig. @ (a) shows 
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FIG. 9. Temperature dependence of Re xi^/c, uji) for cji 
as indicated in Fig. [jj (a) loz — 2. 8 J: result from Eq. (1C). 
(b) Lu 2 = 1.96J and (c) wi = 0.23J: results from Eq. @, 
influence of Auj. For Auj — 0.05 J the points in Fig. scatter 
significantly at intermediate temperatures, the results are re- 
liable for 0.06 J < Alu < 0.15 J. The short-dashed line in (c) 
is the field-theoretical result Eq. (^) deviating significantly for 
ksT < 0.6J. The inset shows the scaling function defined in 
Eq. J2s[ ) for different chain lengths. The frequency at L — 14 
is uji = 0.27J and at L = 12 it is ui\ = 0.31J. 



the rescaled imaginary part of the rescaled field theoreti- 
cal curve, the dashed line in Fig. (a) shows the real part 
at k B T = 0.3 J. Both fit the numerical data better than 
the unsealed result (full lines). Above k B T > 0.3J the 
field-theoretical function Eq. (^|) is not a good fit func- 
tion any more since intermediate frequency values are 
underestimated . 

In the temperature window 0.1J < k B T < 0.3J rele- 
vant for describing the quasi-elastic scattering in CuGeC>3 
the above described fitting procedure is sufficiently accu- 
rate. More precise results for the correlation function 
can be obtained by fitting a spline to the spectrum of the 
real part and simultaneaously fitting its Kramers-Kronig 
transform to the imaginary part. 



VI. APPLICATION TO CUGEO3 

CuGeOs undergoes a spin-Peierls transition at Tsp — 
14.3 K or fceXsp ~ The wave vector of the modu- 

lation in the ordered phase is qo = (it / 'a, 0, 7r/c), where a 
and c are the lattice constant along the crystallographic 
x and z direction, respectively. I set J/k B = 150 K which 
is together with a value of Je/|/ j= 0.24 among those dis- 
cussed as valid for CuGeOs.l 
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FIG. 10. Intensity of quasi-elastic scattering, (a) The fast 
drop I ~ (T— Tsp) -1 above the phase transition at Tsp = 14.3 
K is basically independent of the choice of parameters and 
matches the data from neutron scattering (filled circles, Ref. 

(b) The curve of T/J(qn, 0) obtained with the larger res- 
olution <j x = 0.5 J match the X-ray data (open circles, scaled, 
Ref. \v\ ) better than that obtained with <j„ = 0.05 J. (c) The 
influence of the choice of parameters shown on a logarithmic 
scale. For cr x = 0.5J the curves lie up to an order of mag- 
nitude higher than for o n = 0.05J. The dot-dashed lines 
show the effect of omitting the frequency scaling (gr = 1). 
The dashed curves show the effect of the different functional 
dependence of XoQ^bT/J) in Eq. \n) for a — J2/J = 0. 



A. Quasi-elastic scattering 

Quasi-elastic scattering has been—observed in X-ray 
scattering up to 40 K or k^T w 0.3 J.EZl I can thus use the 
results presented in Sec. [v| to compare with experiment. 
In neutron scattering experiments th e . qu asi-elastic scat- 
tering is observable only up to 16 K.E£rE3 

The frequency dependent scattering rate of inelastic 
neutrons or X-rays can be obtained by convoluting the 
dynamical structure factor (|l]) with a Gaussian of the 
width of the experimental resolution a. 



J(qo,w) 



V2no- J- 



e- 1 ^!^ S*(qoV) duj' (30) 



The energy resolution in neutron scattering is of the 
order of a few meV while X-rays integrate over a much 
larger energy interval, a = oy-p 0.05 J simulates the res- 
olution of diffracted neutronsEi and a = a x ~ 0.5 J is rel- 
evant for X-ray scattering. The X-ray resolution is prob- 
ably even larger, but the interval of —0.5 J < lj < 0.5J 
covers the full width of the relevant magnetic spectrum, 
compare Fig. 0. Also, the fit function Eq. (]|) is not reli- 
able for \ui\ > J. 

Fig. [lO] (a) shows the fast drop of the intensity of the 
quasi-elastic scattering above the phase transition. The 



theoretical intensity (full line) is normalized with Iq — 
-f(qo, 0)| t= 14.4 k- It is basically independent of the choice 
of parameters, all curves shown in Fig. |l^ (c) fall more 
or less on the same line on the linear scale. 

The full curve in Fig. [Io| (a) matches the data from 
neutron scattering (filled circles, Ref. [l5]). Note that 
the critical region has been estimated via the Ginzburg 
criterion^ to be T SP ± 0.4 K. Within this region the 
theoretical divergence of the intensity I — 5(qo,0) ~ 
(T-Tsp)- 1 (dotted line in Fig. |o| (c)) will be suppressed 
by critical fluctuations. A more quantitative analysis 
also must include the limited momentum resolution of 
the neutron data as discussed in Ref. [l4|. 

Fig. |l0| (b) shows a plot of T//(q ,0). As desired, 
the curve obtained with the larger resolution a x = 0.5 J 
match the X-ray data (open circles, scaled, Ref. [l7]) bet- 
ter than that obtained with a n = 0.05 J. 

The influence of the choice of parameters becomes ap- 
parent on the logarithmic scale shown in Fig. [l(] (c). For 
o~ x = 0.5 J the curves lie an order of magnitude higher 
than for a n — 0.05 J at large temperatures. This explains 
why the correlation length can be extracted to higher 
temperatures frpm X-ray measurementai3 then from neu- 
tron scattering.E3 The dot-dashed lines show the effect of 
omitting the frequency scaling (g<r = 1). For a n = 0.05 J 
the intensity is suppressed since the imaginary part with- 
out rescaling has less weight in the low energy regime (see 
Fig. H (a)). For a„ = 0.5 J this effect is cancelled by the 
larger contribution of the real part (Fig. [?] (a)). The 
dashed curves show the effect of the different functional 
dependence of Xo ( k BT/J) in Eq. (§) for a = J 2 /J = OJB 

Overall, the effect of energy resolution accounts for the 
largest effects in the temperature dependence of the in- 
tensity of the quasi-elastic scattering. The effect of the 
frequency scaling gx is a minor correction. Both neu- 
tron and X-ray scattering intensities are satisfactorily de- 
scribed within RPA. 



B. Phonon hardening 

The optical zone-boundary phonons coupling to the 
spin chains in CuGeOa have been, shown to harden 
with the lowering of temperature.Ea The phononic ex- 
citations in RPA are essentially given by the roots of the 
real part o f ^he .1 d enominator of the normal coordinate 



propagator .I13E3IU 



Rex(q,w) = [X) l/ 9 , ! / 1 qfl , ! /,-q D^iq^uj) 



(31) 



The coupling constants and bare phonon frequencies are 
given in Ref. [l2[ The Peierls-active phonon mode cou- 
pling strongest being il2. qo = 2. 13 J, the left hand side 
of Eq. (|3l]) can be approximated by Rex(q, 1-96 J) as 
calculated in Sec. [y]. 

The resulting temperature evolution is up to a pref- 
actor and a constant given by Fig. ^ (b). In contrast to 
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the case of determining the dimer-dimer copelation func- 
tion by the field-theoretical result Eq. @ ,Ej the phonon 
hardens with decreasing temperature while k^T > 0.2J. 
This is reminiscent of the change of sign of the real part 
of the dimer-dimer correlation function discussed in Sec. 
0. The temperature evolution is in qualitative agreement 
with the experimental data from Ref. [l8|, but the phonon 
hardens only by 1% which is less than what is experimen- 
tally observed. The higher Peierls-active phonon modes 
show no spin-phonon coupling induced temperature de- 
pendence both in experiment and theory. 



VII. CONCLUSIONS 

Within the energy resolution limited by finite system 
sizes the dimer-dimer correlation functions obtained by 
exact diagonalization yield satisfactory approximations 
to the thermodynamic limit. This is supported by the 
comparison with the exact results for the XY model. The 
determination of both the real and the imaginary part 
yields an essential improvement of the reliability of the 
results. 

The quantitative applicability of field-theoretical re- 
sults for the dimer-dimer correlation functions is drawn 
into doubt, even at low frequencies and temperatures. 
At intermediate temperatures they can be used as fit 
functions by appropriately rescaling the amplitude and 
frequency dependence. 

The intensity of the quasi-elastic scattering in CuGeC>3 
is satisfactorily reproduced within the RPA approach. 
The impact of the energy resolution of the experimental 
setup on the temperature dependence of the quasi-elastic 
scattering has been shown. The numerical results yield 
the correct qualitative temperature dependence of the 
Peierls-active phonons. 
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